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Abstract 

We compute branching ratios oi B ^ Kir decays in the framework of perturbative QCD 
factorization theorem. Decay amphtudes are classified into the topologies of tree, penguin 
and annihilation, all of which contain both factorizable and nonfactorizable contributions. 
These contributions are expressed as the convolutions of hard h quark decay amplitudes 
with universal meson wave functions. It is shown that (1) matrix elements of penguin 
operators are dynamically enhanced compared to those employed in the factorization 
assumption; (2) annihilation diagrams are not negligible, contrary to common belief; (3) 
annihilation diagrams contribute large strong phases; (4) the uncertainty of current data 
of the ratio R = Br(B° — K^ir^) I'Btl{B^ — K^ir^) and of CP asymmetries is too large 
to give a constraint of ^3. Assuming ^3 = 90" which is extracted from the best fit to the 
data of i?, predictions for the branching ratios of the four B Kir modes are consistent 
with data. 



1 INTRODUCTION 



B fatories at KEK and SLAG are taking data to probe the origin of CP violation. Within 
the Kobayashi-Maskawa (KM) ansatz [|T|, CP violation is organized in the form of a unitarity 
triangle shown in Fig. 1. The angle 0i can be extracted from the CP asymmetry in the B — >■ 
J/ ipKs decays, which arises from the B-B mixing. Due to similar mechanism of CP asymmetry, 
the decays B^ ti^ti^ are appropriate for the extraction of the angle 02- However, these modes 
contain penguin contributions such that the extraction suffers large uncertainty. Additional 
measurements of the decays B^ — > tt^tt" and B^ — * vr^Tr*^ and the use of isospin symmetry may 
resolve the uncertainties [Q]. It has been proposed that the angle 03 can be determined from 
the decays B Kn, rni [|, |, |, Contributions to these modes involve inteference between 
penguin and tree amplitudes, and relevant strong phases have been formulated in terms of 
several independent parameters. We are convinced that one can not make any progress along 
this direction, unless he learns to compute nonleptonic two-body decay amplitudes including 
strong phases. 

The conventional approach to exclusive nonleptonic B meson decays relies on the factoriza- 
tion assumption (FA) [0], in which nonfactorizable and annihilation contributions are neglected 
and final-state- interaction (FSI) effects are assumed to be absent. That is, this approach 
requires simplifying assumptions. Though analyses are easier under this assumption, estima- 
tions of many important ingredients, such as tree and penguin (including electroweak penguin) 
contributions, and strong phases are not reliable. Moreover, it suffers the problem of scale 
dependence ||^. It is also difficult to resolve some controversies in the FA approach, such as the 
branching ratios of the B J/ipK^*^ decays [Q. 

Recently, perturbative QCD (PQCD) factorization theorem for exclusive heavy-meson de- 
cays has been proved, and applied to the semileptonic B D'^*\t{)Ii' decays |Ty, |TT], the 
nonleptonic B — > D^*^7r(p) decays ||, 0] and the penguin-induced radiative B — > K*'j decays 
HTSl , p^ . PQCD is a method to separate hard components from a QCD process, which can 
be treated by perturbation theory. Nonperturbative components are organized in the form of 
hadron wave functions, which can be extracted from experimental data. Here we shall extend 
the PQCD formalism to more challenging charmless decays such as i? ^ Kit, tttt. It will 
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be shown that the difficulties encountered in the FA approach can be resolved in the PQCD 
formalism. 

In this paper we shall evaluate the branching ratios of the following modes, 

Contributions from various topologies, such as tree, penguin and annihilation, including both 
factorizable and nonfactorizable contributions, can all be calculated. That is, FA is in fact 
not necessary. In our approach strong phases arise from non-pinched singularities of quark 
and gluon propagators in nonfactorizable and annihilation diagrams. As explicitly shown in 
Sec. VII, strong phases from the Bander-Silverman-Soni (BSS) mechanism [1^, which is a 



source of strong phases in the FA approach, are of next-to-leading order and negligible. 
As an application, we derive the ratio R and the CP asymmetries defined by 

^ Br(ij° ^ K^TT^) 

Br(5± -> KH^) ' ^ ' 

^ Br(^o ^ K-'K+) - Br(i?g ^ K+n-) 

Br(50 ^ K-TI+) + Br(50 ^ K+n-) ' ^ ' 

^ Bt{B- ^ K^iT-) - Br(i?+ ^ K%+) 
CP Br(fi- ^ i^Oyr-) + Br(fi+ ^ i^07r+) ' ^ 
Vo ^ Br(^g ^ ir°7r°) - Br(i?° ^ ir°7r°) 

Br(50 ^ irovro) Br(50 ^ 7^%°) ' ^ 
^ Br(g- ^ i^-vr") - Br(i?+ ^ X+tt") 

Br(S- ^ ir-vrO) + Br(5+ i^+vrO) ' ^ ^ 

as functions of the unitarity angle 03 using PQCD factorization theorem. In the above expres- 
sions Br(5)^ -> K^-n^) represents the CP average of the branching ratios Br(i?^ K^-k ) and 
Br(i?° K^TT^), and the definition of Br(i?^ — > K^tt^) is similar. It will be shown that the 



uncertainty in the data for R, A^p and A^p |T^, p!7| , 

i? = 0.95 ± 0.30 , A^p = -0.04 ±0.16 , A^p = 0.17 ± 0.24 , (7) 

is still too large to provide useful information of 03. Using the central values of the CLEO data 
for R, we obtain 03 = 90°. 
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An essential difference between the FA and PQCD approaches is that the hard scale at which 
Wilson coefficients are evaluated is chosen arbitrarily as nib or mb/2 in the former, being the 
b quark mass, but dynamically determined in the latter. It has been shown that choosing this 
dynamically determined scale minimizes higher-order corrections to exclusive QCD processes 
||18|| . We observe that this choice leads to an enhancement of penguin contributions by nearly 
50% compared to those in the FA approach. As elaborated in Sec. V, this penguin enhancement 
is crucial for the explanation of the data of all B Kit, vrvr modes using a smaller angle 
03 ~ 90°. Note that an angle 03 larger than 110° must be adopted in order to explain the 
above data in the FA approach ||19|| . 



PQCD factorization theorem for exclusive nonleptonic B meson decays are briefly reviewed 
in Sec. II. The factorization formulas for various B — > Kir decay modes are derived in Sec. III. 
The numerical analysis, including the determination of meson wave functions, is performed in 
Sec. IV. We emphasize the importance of the penguin enhancement in the PQCD approach 
in Sec. V. FSI effects are discussed in Sec. VI. The PQCD approach is compared with other 
approaches in Sec. VII. Section VIII is the conclusion. 

2 FACTORIZATION THEOREM IN BRIEF 

We first sketch the rough idea of PQCD factorization theorem and of its application to two-body 
B meson decays. Take the B ^ tt transition form factor in the fast recoil region of the pion as 
an example. Obviously, this process involves two scales: the b quark mass mf,, which provides 
the large energy release to the fast pion, and the QCD scale Aqcd, which is associated with 
bound-state mesons. Therefore, the B ^ n transition form factor contains both perturbative 
and nonperturbative dynamics. 

In perturbation theory nonperturbative dynamics is reflected by infrared divergences in 
radiative corrections. It has been shown order by order that these infrared divergences can 



be separated and absorbed into a B meson wave function or a pion wave function A 
formal definition of the meson wave functions as matrix elements of nonlocal operators can be 
constructed, which, if evaluated perturbatively, reproduces the infrared divergences. Certainly, 
one can not derive a wave function using a perturbative method, but has to parametrize it as 



a parton model, which describes how a parton (valence quark, if a leading-twist wave function 
is referred) shares meson momentum. The meson wave functions, characterized by Aqcd, niust 
be determined by nonperturbative means, such as lattice gauge theory and QCD sum rules, or 
extracted from experimental data. In the practical calculation below, small parton transverse 
momenta fc^ are included, and the characteristic scale is replaced by 1/6 with b being a variable 
conjugate to k^. 

After absorbing infrared divergences into the meson wave functions, the remaining part of 
radiative corrections is infrared finite. This part, called a hard amplitude, can be evaluated 
perturbatively in terms of Feynman diagrams for decays of an on-shell b quark. Note that 
the b quark carries various momenta, whose distribution is described by the parton model 
introduced above. The hard decay amplitude is characterized by the virtuality t of involved 
internal particles, which is of order m^. The B ^ tt transition form factor is then expressed 
as the convolution of three factors: the B meson and pion wave functions, and the hard b 
quark decay amplitude. This is so called factorization theorem. Note that the separation 
of nonperturbative and perturbative dynamics is quite arbitrary. This arbitrariness implies 
that a renormalization-group (RG) improvement of the factorization formula for the B ^ % 
transition form factor can be implemented. The RG evolution from the all-order summation of 
large logarithmic corrections to the above convolution factors will be made explicit below. 

A salient feature of PQCD factorization theorem is the universality of nonperturbative wave 
functions. Briefly speaking, the infrared divergences associated with the B meson are process- 
independent, and the formal definition of the B meson wave function in terms of matrix elements 
of nonlocal operators is universal for all B meson decay modes. It is not difficult to understand 
this universality: infrared divergences correspond to long-distance effects, while the hard b quark 
decay occurs in a very short space-time. It is natural that these two dramastically different 
subprocesses decouple from each other. That is, the long-distance dynamics is insensitive to 
specific decays of the b quark with large energy release. Because of universality, a B meson 
wave function extracted from some decay modes can be employed to make predictions for 
other modes. This is the reason PQCD factorization theorem possesses a predictive power. 
We emphasize that PQCD is a theory, instead of a model, since higher-order and higher-twist 
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contributions can be included systematically. The model independence of PQCD predictions 
can be achieved, once wave functions are determined precisely. 

PQCD factorization theorem for nonleptonic B meson decays, such as B ^ K{ti)ti and 
is more complicated. These decays involve three scales: the W boson mass 
Mw, at which the matching conditions of the effective weak Hamiltonian to the full Hamiltonian 
are defined, the typical scale t (of order rrih) of a hard amplitude, which reflects the specific 
dynamics of a decay mode, and the factorization scale 1/6 (of order Aqcd) introduced above. 
The dynamics below 1/6 is regarded as being completely nonperturbative, and parametrized 
into meson wave functions 0(x, 6), x being the momentum fraction. Above the factorization 
scale, radiative corrections to on-shell 6 quark decays can be computed perturbatively. This part 
contains two characteristic scales, M\y and t, differing from the case of the B it transition 
from factor, and further factorization is necessary. 

Radiative corrections produce two types of large logarithms: ln(Mvi//t) and ln(t6). The 
former are summed by RG equations to give the evolution from Mw down to t described by the 
Wilson coefficients C(t), while the latter are summed to give the evolution from t to 1/6. The 
matching between the full Hamiltonian and the effective Hamiltonian in the above three-scale 
factorization theorem is similar to that in the standard effective field theory. The difference is 
that diagrams in the full theory contain not only W boson emissions but hard gluon emissions 
from spectator quarks One can show that the effective operators, in the presence of the hard 
gluons from spectators, still form a complete basis, and that the Wilson coefficients derived in 
the three-scale factorization theorem are the same as those derived in the standard effective 
theory. 

There also exist double logarithms ln^(P6) from the overlap of two types of infrared di- 



vergences, coUinear and soft, in radiative corrections to meson wave functions pO|, where P 



denotes the dominant light-cone component of a meson momentum. The resummation pO| , pi 



of these double logarithms leads to a Sudakov form factor exp[— s(P, 6)], which suppresses the 
long-distance contributions from the large 6 region, and vanishes as 6 = 1/Aqcd- This factor 
guarantees the applicability of PQCD to exclusive decays around the energy scale of the 6 
quark mass |jTO|. For a detailed derivation of the relevant Sudakov form factors, we refer the 
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readers to []TU|, |TT|. With all the large logarithms organized, the remaining finite contributions 
are absorbed into the hard b quark decay amplitude H(t). 

A three-scale factorization formula for exclusive nonleptonic B meson decays possesses the 
typical expression, 



C(t) (g) H{t) ® (t>{x, h) (g) exp 



-.(P,6)-2 f %ia,m 

Jl/b LL 



(8) 



1/6 /i 

where the exponential involving the quark anomalous dimension 7 = —Us/t^ describes the 
evolution from t to 1/6 mentioned above. Note that Eq. (§) is a convolution relation, with 
internal parton kinematics x and h integrated out. The hard scale t, related to the virtuality of 
internal particles in hard amplitudes, depends on x and h. All the convolution factors, except for 
the wave functions 6), are calculable in perturbation theory. The wave functions, though 
not calculable, are universal. If choosing t as the h quark mass m;,, the Wilson coefficient 
C(mfe) is a constant, and Eq. reduces to the simple product of the Wilson coefficient and 
a hadronic matrix element. However, the analysis of next-to-leading-order corrections to the 



pion form factor |T8| has suggested that t should be chosen as the virtuality of internal particles 
in order to minimize higher-order corrections to the hard amplitudes. 

Because of the universality of nonperturbative wave functions stated above, the strategy 
of PQCD factorization theorems is as follows: evaluate all perturbative factors for some decay 
modes, and adjust the wave functions such that predictions from the corresponding factorization 
formulas match experimental data. At this stage, the nonperturbative wave functions are 
determined up to the twist, at which the factorization is constructed. Then evaluate all the 
perturbative factors for another decay mode. Input the extracted wave functions into the 
factorization formulas of the same twist, and make predictions. With this strategy, PQCD 
factorization theorems are model-independent and possess a predictive power. In Sec. IV we 
shall make explicit the determination of the B meson, kaon, and pion wave functions from 
currently available data and phenomenological arguments. 
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3 B ^ Kn AMPLITUDES 

The effective Hamiltonian for tlie flavor-clianging b ^ s transition is given by 

10 

Ci(/i)oi^)(/i) + C2(/i)o(^V) + ECi(/i)o.(/i) 



V ^ q=u,c 



1=3 



(9) 



witli tlie Cabibbo-Kobayaslii-Maskawa (CKM) matrix elements Vg = V*Vgb and tlie operators 





= {Siqj)v-A{qjbi)v-A , 






V-A{qjbj)v-A , 


03 = 


{sibi)v-AY{qjqj)v-A , 
q 


O4 = (s, 




-A'Y{qjqi)v-A , 

q 


05 = 


{sibi)v-AY{qjqj)v+A , 

q 


Oe = (s, 




-AY{qjqi)v+A , 

q 


Or = 


3 

-is^bi)v^A J2 (^qimjh+A 

^ q 


, 08 = 




'^j)v-Aj2'^q(^j9i)v+A 

q 


09 = 


3 

-{S^bi)v-A J2 ^qimjh-A 

^ q 


, Oio = 




ibj)v-AYeg{qjqi)v-A 
q 



(10) 

i, j being tlie color indices. Using the unitarity condition, the CKM matrix elements for the 
penguin operators O3-O10 can also be expressed as + K = —Vt- We define the angle 03 via 



Vub = |Kfe|exp(-i03) . 
Here we adopt the Wolfstein parametrization for the CKM matrix upto 0{X^) 



(11) 



Vcd Vcs Vcb 

Wtd Vts VtJ 



2d 
2 



-A 1 

\AX^{1- p-ir]) -AA2 



1 



(12) 



A recent analysis of quark-mixing matrix yields p3| 



A 
A 
Rb 



0.2196 ±0.0023 
0.819 ±0.035 , 



p^ + T]^ = 0.41 ±0.07 . 



For the K^tt^ decays, the operators Oi^l contribute via an annihilation topology, 

and Oi^2 do not contribute at leading order of a^. The absorptive part of the charm quark loop 
integral computed by BSS is thus of higher order. 03,4,5,6 contribute via tree and annihilation 
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(13) 



topologies, and the tree topology involves the B ^ ir form factor. 03,5 gives both factorizable 
and nonfactorizable (color-suppressed) contributions, while 04,6 gives only factorizable ones 
because of the color flow. The contributions from 07,8,9,10 are the same as 03,4,5,6 except for an 
additional factor (3/2) with the light quark q — d in the tree topology and with q — u in the 
annihilation topology. For the B^ K^n^ decays, the operators Oi^2 contribute via a tree 
topology, and Oi^2 do not contribute at leading order of a^. The penguin operators contribute 
in the same way as in the B"^ ^ K^tt"^ decays but with the light quark q — u in the tree 
topology and with q — d in the annihilation topology. The lowest-order hard b quark decay 
amphtudes are summarized in Fig. 2 for B^ — > K^n^ decays and in Fig. 3 for B^ — > K^n^ 
decays. 

For the B^ K^n^ decays, the operators 0^1 contribute via tree and annihilation topolo- 
gies, where the tree topology involves both the S — > tt and B ^ K form factors. The penguin 
operators also contribute via tree and annihilation topologies with the light quark q — uin both 
topologies. While the tree topology involves only the S — > tt form factor. For the B^ — > ir°7r° 
decays, the operators Oi^^ contribute via the tree topology, which involves only the B ^ K 
form factor. The penguin operators contribute in a similar way but with the light quark q = d 
in both the tree and annihilation topologies. Their lowest-order diagrams for the hard h quark 
decay amplitudes are basically similar to those in Figs. 2 and 3. 

The momenta of the B and K mesons in light-cone coordinates are written as Pi = 
(Mb/v/2)(1, 1,0t) and P2 = (Mb/V2)(1, 0, Oy), respectively The B meson is at rest with 
the above parametrization of momenta. We define the momenta of light valence quark in the B 
meson as ki, where ki has a plus component kf, giving the momentum fraction xi — kf/Pi, 
and small transverse components ki^. The light valence quark and the s quark in the kaon 
carry the longitudinal momenta X2P2 and (1 — X2)P2, and small transverse momenta and 
— k2T, respectively. The pion momentum is then P3 = Pi — P2, whose nonvanishing component 
is only P3". The two light valence quarks in the pion carry the longitudinal momenta X3P3 and 
(1 — a;3)P3, and small transverse momenta ii^T and —\s.3t, respectively. The kinematic variables 
associated with each meson are indicated in Fig. 4. 

The Sudakov resummations of the large logarithmic corrections to the B, K and tt meson 
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wave functions (pB, <Pk and (p-,, lead to the exponentials exp(— Sb), exp(— S'a') and exp(— ^tt), 
respectively, with the exponents |jlO|, ^ 

SBit) = s(xiP+,6i) + 2 r ^7(«,(/i)), 

Jl/bi jl 

SK{t) = s{x2P^M) + s{{l-X2)P^M) + 'i f — 7(as(/i)), 

"'1/62 

S^{t) = s{x3P,-,h) + s{{l-X3)P,-,h)+2 f (14) 

Ji/b-i fi 

The variables bi, 62, and 63 conjugate to the parton transverse momentum kiT, k2T, and ksT 
represent the transverse extents of the B, K, and vr meson, respectively. The exponent s is 



written as 20, 21 



Ji/b fi 



(15) 



where the anomalous dimensions A to two loops and B to one loop are 



A 



a. 



Cf— + 



67 10 2 /e^- 



3 TT I 2 



(16) 



with Cp = 4/3 a color factor, f = 4: the active flavor number, and 7^; the Euler constant. The 
one-loop expression of the running coupling constant. 



"^^^^-/?oln(/iVA^) ' 
is substituted into Eq. (p!5|) with the coefficient Pq = (33 — 2/)/3. 
The decay rates of B^ K^tt^ have the expressions 

„ GlMl 



;i7) 



1287r 



:i8) 



The decay amplitudes A'^ and ^ corresponding to B^ K^n^ and i? K^n , respectively, 
are written as 



A' 
A- 



fKv;F[ + v:mi + fBv:F^ + v:mi - fBv:F, - v:Ma 

fKVtF[ + VtMl + fBVtF^ + VtMl - fsVuFa - V^Ma , 



(19) 
(20) 
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with the B meson (kaon) decay constant fB(K)- The notations F represent factorizable con- 
tributions (form factors), and M. represent nonfactorizable (color-suppressed) contributions. 
The indices a, e and P denote the annihilation, tree and penguin topologies, respectively. F^, 
associated with the time-like K ^ it form factor, and Aia are from the operators O^j^^. The 
decay rates of i?^ K^ti^ have the similar expressions with amplitudes 

A = fKV:Ff + V;M^ + IbV^F,'' + V;M^ - fKV:F,-V:Me , (21) 
A = fKVtF[ + VtM^ + fBVtF[ + VtM^-fKVuF,-V^Me, (22) 

for i?J] — * K^TT^ and 5^ —>■ K^n^, respectively. The notations are similar to those in Eqs. (p!9|) 
and (ppp. Fg, associated with the i? — > vr form factor, and A^g are from the operators Oi^2- 
The decay amplitudes for — > K^tt^ are given by 

72^'+ = fKv;F[ + v;m!: + fBv;F[ + v;Mj: - fBv:F,-v:Ma 

-fKv:F, - v:Me - uv:f,k - v:m.k , (23) 

v^^'- = fKVtF!^VtM'',^jBVtFl^VtMl-fBVuFa-VuMa 

-fKVuFe - VuMe - f.VuF.K ' VuM^K , (24) 

which correspond to B^ — > K^tt^ and B^ — > K^n^, respectively. The factorizable contribution 
Fex is associated with the B K form factor, and AieK is the corresponding nonfactorizable 
contribution. Similarly, the decay rates of B^ K^n^ are obtained from the amplitudes 

V2A' = fKV;Ff + v:M^ + fBV:F[ + v;M^-f.V:F,K-V:MeK, (25) 
V2A' = fKVtF[ + VtM^ + fBVtF^ + VtM^ - UVuF.K -VuMeK , (26) 

for B^ A'°7r0 and B^ K^7r°, respectively. 

Basically, one needs to derive the factorizaton formulas only for the tree and annihilation 
topologies. Wilson coefficients corresponding to different operators are then inserted into the 
factorization formulas. The form factors are written as 

pP - pP + pP 

/•l roc 

= IQ-nCpM"^ / dxidxs / bidbib^dbscpBixi^bi) 
Jo Jo 
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X{[(1 +X3)0,r(a;3) - 2X3)(f)'^{x3)] Ee4{t'e^)he{Xi, X3,bi,b3) 

+2r^<f)'^{x3)E,4{tf^)K{x3,Xuh,h)} , (27) 

= 32nCpM'^ / dxidxs / bidbib3db3(l)B{xi,bi) 
Jo Jo 

y<rK{[(pn{x3) + r^(2 + X3)0^(X3)] Ee6(4^^)/ie(Xi, X3, 61, 63) 

+ [x^M^s) + 2r.(l - a;i)0;(2;3)] ^e6(4'^)/ie(x3, ^i, 63, ^l)} , (28) 



pP - pP > pP 

1 i"00 



pi poo 

= IQnCpM^ / dx2dx3 / 62'^^2&3'^&3 

x{[-X3(I)k{x2)(PA^3) - 2r^rK{l + X3)(f)'j^{x2)(p'^ix3)] Ea4it'^^^)ha{x2, X3, 63, ^s) 

+ [x2(pK{x2)(f)7T{x3) + 2r^rx(l + X2)(P'k{x2)(P'^{x3)] EaA{t^a^)ha{x3, X2, 63, ^2)} , (29) 

F[q = 327rC^M| / dx2dx3 / 62t^fe2&3t^&3 

JO JO 

x{[r^X30x(a;2)0Ua;3) + 2rK(p'K{x2)(pn{x3)] EaQ{t^a^)ha{x2,X3,b2M) 

+ [2r^0A'(x2)0Ua;3) + rKX2<p'K{x2)4>A^:i)\ Ea&{tf^)ha{x3, X2, 63, ^2)} , (30) 

with the evolution factors 

E,i{t) = asit)aiit)exp[-SBit)-S^it)] , (31) 
Eaiit) = as{t)ai{t)exp[-SK{t)-S^{t)]. (32) 

The expression of F^ (Fa) for the O12 contributions is the same as F^ (-^^) but with the Wilson 
coefficient ai(te) (a.i(ia))- The factorization formula of F^k is written as 

pi poo 

FeK = IQirCpM"^ / dxidx2 / bidbib2db2(t)Bixi,bi) 
Jo Jo 

X{[(1 + X2)(j)K{x2) + rx(l - 2X2)0x(a;2)] -E2(tix)/ie(Xl,X2,&l,&2) 
+2r;^0^(x2)F2(tf/J)/ie(x2, Xi, 62, &l)} , (33) 

with the evolution factor 

E2{t) = a,{t)a2{t) eM-Ssit) - SK{t)] . (34) 

The hard functions /I's in Eqs (^7|)-(^) and in Eq. (p3D, are given by 

/ie(a;i,X3,6i,63) = KQ{^/x^MBbi) 

12 



X [e{hi - h)Ko (v^Mb6i) Jo (v^M^fes) 
+e{h-h)Ko{V^3MBh)Io{V^3MBh)] , 



(35) 



X 



(1) 



e{b2 - h)K' (v^Msfos) Jo (v^Mb63 
+^(63 - b2)Hi,''> (V^Msfes) Jo (V^Mb^s)" 



(36) 



The derivation of h, from the Fourier transformation of the lowest-order H, is the same as that 
for the B Dtt decays [^], but with a vanishing D meson mass. The hard scales t are chosen 
as the maxima of the virtualities of internal particles involved in b quark decay amplitudes, 
including l/bf. 



''e 


= max(i 


v^Ms, 1/61 


A/bs), 


+(2) 
''e 


= max(i 


/^Mb, 1/61 


, I/&3) , 




= max(i 


/^Ms, 1/62 


A/bs), 




= max(i 


/^Mb, 1/62 


, l/bs) , 


''eK 


= max(i 


/^Mb, l/bi 


,l/&2), 


''eK 


= max(i 


^iMb, 1/61 


, I/&2) , 



(37) 



which decrease higher-order corrections. The Sudakov factor in Eq. (14) suppresses long- 
distance contributions from the large b region, and improves the applicability of PQCD to 
B meson decays. 

For the nonfactorizable amplitudes, the factorization formulas involve the kinematic vari- 
ables of all the three mesons, and the Sudakov exponent is given hj S = Sb + Sk + S-,,. The 
integration over 63 can be performed trivially, leading to &3 = &i, 63 = 62, or 62 = ^i- Their 
expressions are 



Mi 



e4 



?,2t:Cf^2N^MI ^ [dx] bidbib2db2^Bixi, 6i)0i^(x2) 

x{[(2;i - X2 - X3)0^(x3) + r^X3(j)'^{x3)] E'^^{t^J'')h^J\xi, X2, X3,bi,b2,bi) 

('2) ('2) 

+ [(1 -Xi -X2)0,r(a;3) - r^X3(f)'Ax3)] Ei,4^{ty)hy {Xi, X2, X3,bi,b2,bi)} 



(38) 
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>^rK{[{xi - X2)(l)^{x2) + r^{xi -X2- x^)(t)'^{x^)\ E'^f^{t^p)h^^\xi, X2, X3, 61, 62, &i) 
+ [(1 -xi- X2)(t)^{x^) + r^(l -X1-X2 + x^)(t)'^{x^)] 

xE'^Sf)hf{x^,X2,x^MMM)} , (39) 

M^^ = 2,2tiCfsJ2N,MI j^[dx\ hidhih2dh2<j)B{xiM) 

x{[xz(t>K{x2)(t>T:{x^) - r^rxixi - X2 - X3)(f)'i^{x2)(l)'^{x3)] 

xE'^^{tf)hf{x,,X2,x^MMM) 

- [{Xl + X2)(t)K{x2)(t>-K{x3) + Vk(2 + Xi + X2 + X3)0';^(x2)0^(x3)] 

x£;;4(if)/lf(xi,X2,X3, 61,62, &2)} , (40) 

= ^2tiCfJ2N,MI f\dx] H hidhih2dh2(l)B{xiM) 

x{[-r^X30/f (X2)0UX3) - rK{xi - X2)^'k{x2)^'k{x2)] E'^Q{tf^)hf^ {Xi, X2, X3, bi, 62, 62) 

- ['^7r(2 - X3)0;f (X2)0^(X3) - T/f (2 - Xi - ^2)0';^ (x2)07r(a:^3)] 

x^;6(tf )/if (a:i,a;2,a;3,6i,62,62)} , (41) 
with the number of colors N^. = 3, the definition [dx] = dxidx2dx3, and the evolution factors 

E'Jt) = a,(t)<(t)exp[-^(t)|,3=,J, (42) 
E'^,{t) - a,{t)a\{t)exp[-S{t)\,,=b,]. (43) 

The expression of A4e {A4a) is the same as i^ad) but with the Wilson coefficient a[{td) 
{a'i{tf)). The nonfactorizable amplitude A4eK is written as 

I rl roc 

MeK = 327rCF^2N,MljJdx]j^ hidhihdh(t)B{xiM)(t)Ax^) 

X{[(xi -X2 -X2)4)k{x2) +rKX2(p'Kix2)]E'2{t''^^)h'^j\xi,X3,X2,bi,bi,b3) 

+ [(1 -Xi- X3)(I)k{x2) - rKX2(p'K{x2)] E2{t^^}^)h''^\xi, X3, X2MM,h)} , (44) 

with the evolution factor 

&2=&1 ] ■ (45) 
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The functions h^^\ J = 1 and 2, appearing in Eqs. (|5BD-(^TD and in Eq. (PD, are written as 



[9ih - b2)Ko (DMsh) Jo {DMBb2) 

+e{h - h)Ko {DMBh2) lo {DMBhi)] 
xKo{DjMBb2) , 



in r 
'2 



mMBb2 



for D]>0, 
for < 



(46) 



9{h - h2)H^^^ (FMBh) Jo (FMB62) 
+^(62 - fei)i^i'^ iFMBb2) Jo {FMbW) 
xKoiFjMBh) , 

X 



for F/ > 



F/|Ms6i) 



for F/ < 



(47) 



with the variables 



?2 



Fi = {Xi - X2)X3 , 
-(1 - Xi - X2)X3 , 

Xi + X2 + (1 - Xi - X2)X3 . 

For details of the derivation of h^^\ refer to ||2^. The hard scales t^^^ are chosen as 



^2 



F^ 



(48) 





= m8ix{DMB, 


1Mb, 1/61 


I/&2) , 


^d 


= m8ix{DMB, 


v/|Di|MB,l/6i 


I/&2) , 




= max(FMB, ^ 


/|F2|Mb,1/6i, 


I/&2) , 


.(2) 


= max(FMB, a 


J\Fi\MB,l/h, 


I/&2) , 


'^dK 


= max(DMB, 


^\DI\MbMW 


I/&3) , 


^dK 


= max{DMB, 


J\Dl\MBMh 


I/&3) . 



(49) 



In the above expressions the Wilson coefficients are defined by 



ai 
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a, = C, + - 



t ^2 

ar. - 



a, 

C'3 3 / Cg 



If 3 



= ^ (Cs + ^e.C^) . (50) 

Both QCD and electroweak penguin contributions have been included as shown in Eq. (|50|). 
The factors and Vk, 

mo^ Ml 
Mb rriu + rud 

moK Ml 
Mb rus + nid 

with m^, mrf, and being the masses of the m, d and s quarks, respectively, are associated with 
the normalizations of the pseudoscalar wave functions (j)' . The pseudovector and pseudoscalar 
pion wave functions and 0^ are defined by 



[X = / — — e 



o(0|%^)7"75«(0)|7r), (53) 



27r 2 

X) = / — e -3"^^i(0|%+)75M(0)|7r), (54) 



3 

satisfying the normalization 

- dx^^{x) = j dx(t)'^{x) = 7-^= , (55) 
Jo l\jzNc 

with the pion decay constant j^^. The kaon wave functions (\)k and 0^ possess similar definitions 
and normalizations with the d quark field, m^, M^r and replaced by the s quark field, m^, 
Mk and /i^, respectively. 

Note that we have included the intrinsic h dependence for the heavy meson wave function 
^B but not for the light meson wave functions ^.^ and ^k- It has been shown that the intrinsic h 
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dependence of the light meson wave functions, resulting in only 4% reduction of the predictions 



for the B ^ n form factor, is not important [|T0|- It is reasonable to assume that the intrinsic 
b dependence of the kaon wave function, which is still unknown, is not essential either. As 
the transverse extent b approaches zero, the B meson wave function (f)B{x,b) reduces to the 
standard parton model 4>b{x), i.e., 0_b(x) = (pB{x,b = 0), which satisfies the normalization 



fB 



'^^^^^^^ = 2vm • ^^^^ 

We do not distinguish the pseudovector and pseudoscalar components of the B meson wave 
functions under the heavy quark approximation. They have roughly the same normalizations 
because of MB/{mb + nid) ~ 1. 



4 NUMERICAL ANALYSIS 



In the factorization formulas derived in Sec. IV, the Wilson coefficients evolve with the hard 
scale t that depends on the internal kinematic variables Xj and 6j. Wilson coefficients at a scale 
fi < My/ are related to the corresponding ones at yU = M\y through usual RG equations. In 
our analysis we adopt the leading-order expressions for the Wilson coefficients with QCD and 
electroweak penguins included. 



exp 



■C{M, 



w 



(57) 



where the leading-order anomalous dimension matrices 7*^°^ are referred to The matching 
conditions at yU = Mw and the choices of the relevant parameters are given in Appendix 
A. 

Since the typical scale t of a hard amplitude is smaller than the b quark mass rrib, we further 
evolve the Wilson coefficients from ^ = rrii, down to n = t using the RG equation. 



Q^s (A'-) ^ (0)T I O^emifJ') ^ (o)r 



(58) 



ATT " 4:71 

where the anomalous dimensions 7^°j for / = 4 are referred to [Q. The solution to Eq. (^) 
and the values of the Wilson coefficients Ci{mb) are also listed in the Appendix A. For the scale 
t below the c quark mass nic = 1.5 GeV, we still employ the evolution function with / = 4, 
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instead of with / = 3, for simplicity, since the matching at rric is less essential. Therefore, we 
set / = 4 in the RG evolution between t and 1/6 governed by the quark anomalous dimension 

7- 

For the B meson wave function, we adopt the model 



')b{x, b) = Nbx'^{1 — xY exp 







\ ^bO 






1 2 



(59) 

with the shape parameter oob = 0.4 GeV 1^^. The normalization constant Nb, which is related 
to the decay constant J'b, will be determined below. As to the pion wave functions, we employ 
the models 

Mx) = ^|=Ax(l-x)[l + c^(5(l-2x)2-l)] , (60) 
= ^/.a:(l-a:)[l + c;(5(l-2a:)2-l)], (61) 

with the shape parameters c^^ and c^. The kaon wave functions are chosen as 

Mx) = -|=/^x(l -x)[l + 0.51(1- 2a;) + 0.3(5(1 - 2x)2 - 1)] , (62) 
V 2iVc 

M^) = ^/i.x(l-x)[l + 4(5(1 -2x)2-l)]. (63) 



(pK is derived from QCD sum rules ||2^, where the second term 1 — 2x, rendering (f)^ a bit 
asymmetric, corresponds to SU{3) symmetry breaking effect. The decay constant fx is set to 
160 MeV (in the convention = 130 MeV). Since predictions for the B Ktt decays are 
insensitive to the kaon wave functions, we simply adopt the result of QCD sum rules. For the 
same reason, we assume that 0'^ and 0^ possess the same functional form and that the shape 
parameter c'^ of the term 5(1 — 2x)^ — 1 in 0'^ is equal to c^. 

We propose to determine from the branching ratios of the B Dn decays: 



- MBf^Tm^ ' ^^^^ 



because this quantity is insensitive to mo7r and 0^. In order to render PQCD predictions reach 
the central value of the data of Rd = 1-61 [^, a large c^r = 0.8, which enhances nonfactorizable 
contributions to the B^ D^n^ decay, is preferred. On the other hand, the data of the 
B pn decays also imply a large Cj^. To further enhance nonfactorizable contributions relative 
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to factorizable ones, the B meson wave function with 6 — i> has been assumed 
as shown in Eq. (|59|) . This behavior, different from that of the model 0b ^Jx as x — 
proposed in decreases factorizable contributions. Note that nonfactorizable contributions 
are insensitive to the variation of the B meson wave funciton. The details for the above 
numerical study will be published elsewhere. The extracted pion wave function 0^ with c^r = 0.8 
is close to the Chernyak-Zhitnitsky model with c^r = 1.0 It differs from the asymptotic 
model with c^r = 0, which has been extracted from the data of the pion transition form factor 



involved in the process 7r7* — > 7 We shall argue that 7r7* ^ 7 is a special process, whose 



infrared divergences differ from those of processes containing final-state hadrons [Q. Hence, 
there is no contradiction between determined from the B — > Dn decays and from the pion 
transition form factor. 

We then extract from the data of the pion form factor, whose factorization formula is 
written as 

Ftt{Q'^) = IQttCfQ'^ / dxidx2 / bidbib2db2as{t) expl—Sj^jrit)] 

Jo Jo 

XX2 [(p7r{Xi)4>n{x2) + 2r^0'(xi)0^(x2)] h{xi, X2, &1, ^2) , (65) 



with 



SnAt) = s{XiP+, 61) + S((l - Xi)P+ , 61) + s{x2P-2, 62) + s((l - X2)P-2, h 



•2 



+2/* %(asm)) + 2f ^7(«.(/^)), (66) 

Jl/bi /i Jl/b2 fJ, 

h{xi,X2,bi,b2) = Ko{./x^Qbi) 

X [^(61 - 62)^^0 iV^Qbi) lo (v^Q&2) 

+e{b2 - b,)Ko {V^2Qb2) lo (v^Q&i)] , (67) 

t = max(-yxiX2<5, 1/&1, 1/&2) • (68) 

The momentum transfer is defined by = 2Pt^i ■ Pt^2i P-kI and Pt^2 being the momenta of the 
intial and final pions, respectively. Useful references for the derivation of the above expression 
are [||, |3|, |35|. The data are Q^F^{Q^) ~ 0.4 ± 0.2 GeV^ for > A GeV^ |6|, g^. Adopting 
the quark masses rriu = 4.5 MeV and rriii = LSnid = 8.1 MeV, which lead to itlq.^ = 1.53 GeV, 
we find that the choice = gives the pion form factor Q'^F^^IQ'^) ~ 0.4 GeV^ for = 6.3 
GeV^. Hence, we choose = = as stated before. 
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With the pion wave functions fixed in the above procedures, we determine the B meson 
decay constant J'b (or A''^ in Eq. (^)) from the data of the B^^ tt^tt^ decay |17|, 



Br(5° ^ TT^TT^) = (4.3+}:^ ± 0.5) x 10"^ , (69) 

where Bt{B^ tt^tt^) represents the CP average of Br(i?^ tt+tt^) and Bt{B'^ 7r^7r+). 
We employ Gp = 1.16639 x 10"^ GeV^^ the Wolfstein parameters A = 0.2196, A = 0.819, 
and Rf, = 0.38, the masses Mb = 5.28 GeV, and the B^ {B~) meson hfetime r^o = 1.55 ps 
{tb- = 1.65 ps) p3[. For the factorization formulas of the B ^ titt decays, refer to For the 
angle 03 = 90°, the result is fs = 190 MeV, which corresponds to Bi{B^ Tr^'^Tr'^) = 4.5 x 10^® 
and the B ^ n transition form factor 

= 0) = 0.25 , (70) 

which is also reasonable. Here q stands for the momentum carried away by the external W- 
emission. The value of is close to that adopted in the PQCD studies of the B Dn and 
B K*'~f decays [|13], and consistent with those from lattice calculations and from 



QCD sum rules |^0| in the literature. The motivation to choose 03 = 90° will be explained 
later. 

We emphasize that the decay constant Jb can not be determined unambigiously in the 
current analysis. The above value Jb = 190 MeV corresponds to the shape parameter ujb = 0.4 
GeV. Changing ujb, different J'b will be obtained when fitting PQCD predictions to the data 
in Eq. (^). However, if more data, such as the CP asymmetry in the tt^tt^ decays, are 

available, both ujb and can be uniquely determined. The reason is that tree and penguin 
contributions depend on ujb and Jb simultaneously, while annihilation contributions, the most 
important source of strong phases as shown below, depend only on J'b- Because the branching 
ratio, mainly determined by tree and penguin contributions, and the CP asymmetry, related to 
strong phases of annihilation contributions, vary with the B meson wave function in a different 
way, their data can fix Ub and Jb uniquely. 

Note that the above parameters are obtained by fitting predictions to the central values 
of the available data. If taking into account the uncertainty of the data, the allowed range of 
the parameters is in fact huge. For example, any value of the shape parameter in the pion 
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wave funciton 0,^ between 0.4 and 1.0 is acceptable for the data of Rd- The shape parameter 
in the pseudoscalar kaon wave function 0^ can differ from in 0^. In this work we do 
not intend to determine the range of parameters, but adopt representative parameters to make 
predictions for the B Ktt decays, and examine whether the predictions are consistent with 
the data. For a summary of the parameters we have adopted in the numerical analysis, refer 
to Appendix B. 

With all the meson wave functions fixed, we predict the branching ratios and the CP asym- 
metries of the B — > Kit decays. The s quark mass is set to = 100 MeV, which corresponds 
to rrioK = 2.22 GeV. We derive the branching ratios of the four B Kit modes in Eq. (|^) for 
different 03, which are shown in Fig. 5. The branching ratios of the K^n^ and K^tt^ modes 
increase with 03, while those of the K^ir^ and K^tt^ modes are insensitive to the variation 
of 03. The increase with 03 is mainly a consequence of the inteference between the penguin 
contribution and the tree contribution Fg. Predictions for the ratio R in Eq. (|^) and the CP 
asymmetries Acp in Eqs. (|])-(|]) for different 03 are displayed in Fig. 6 and Fig. 7, respectively. 
The prediction of R increases from 0.7 to 1.2 when 03 moves from to 180°. Unfortunately, 
the large uncertainty of the current data does not give a constraint of 03. Comparing with the 
central value of the CLEO data of R in Eq. (|^, we extract 03 = 90°. The data of Aqp have also 
large uncertainties, and do not constraint 03 either. Our analysis shows that the magnitude of 
Aqp and A^p is negligible, smaller than 3%, while the magnitude of A^^p and A^p can reach 
20%. 

Our predictions for the branching ratio of each mode corresponding to 03 = 90°, 



Br(5+ - 




= 20.22 X 10 


Bt{B- - 


K\-] 


= 19.79 X 10 


MB', - 


K+n-] 


= 22.74 X 10 


MB', - 


K-TT+] 


= 15.50 X 10 


Br(fi+ - 


K+u^] 


= 11.40 X 10 


Bi{B- - 


K-n'] 


= 7.89 X 10" 


MB', - 




= 8.81 X 10~^ 
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Bt(B 







9.25 X 10" 



(71) 



are consistent with the CLEO data 16 



Br{B^ - 




= (18.2 


tt'o ± 1.6) X 


10- 


Br(5° - 


^ K^n^] 


= (17.2 


tli ± 1.2) X 


10" 


Br{B^ - 


K^n"^] 


= (11.6 


t^rtli) X 10 


-6 


Br(5° - 




= (14.6^ 


1;^^ X 10- 


-6 



(72) 



To derive Eq. (|fl|), we have distinguished the pion masses M^± = 139 MeV and M^-o = 135 
MeV. 

The PQCD resuhs of each form factor and nonfactorizable amphtude involved in the — »• 
K~^7T~ decay are hsted in Table I. It indicates that nonfactorizable contributions are only few 
percents of factorizable ones. This is the reason FA works well for most two-body B meson 
decay modes. However, there are exceptions. For modes whose factorizable contributions 
are proportional to the small Wilson coefficient 02, such as B ^ J/iIjK^*\ nonfactorizable 
contributions become important. Similarly, the term F^k, proportional to 02, is small. Hence, 
the branching ratios of the Ktt'^ modes are about half of those of the Kn^ modes. Table I also 
indicates that the factorizable annhilation diagrams contribute dominant strong phases. The 
reason has been discussed in [^. If expressing the amplitude of the B^ K'^tt~ decay as 



A ~ V;Pe'^^ - V*T, 



(73) 



with the penguin contribution P = IfxFf + fsF^l and the tree contribution T = If^F^l, the 
strong phase (5p is as large as 

6p = 144° . (74) 
This result is consistent with the conclusion drawn from a global fit to data of two-body 



charmless B meson decays , where the strong phase was introduced as a free parameter. 

To test the sensitivity of our predictions to different choices of model wave functions and 
parameters, we have varied the shape parameter ub for the B meson wave funciton from 0.3 to 
0.5, the shape parameter c'^ for the kaon wave funciton from to 0.8, the masses mQK{-K) from 
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1.3 GeV to 2.7 GeV, the forms of the meson wave functions, such as 



= iV^^^^x(l - exp 













' 2 



(75) 



4>f\x) = + 0.3(5(1 -2x)2-l)], (76) 

for the B meson and the kaon, and the asymptotic model 

0f (^) = - ^) ' (77) 

for the pion, and the Wolfstein parameters A from 0.21 to 0.22. It is found that our predictions 
for R change by less than 5%, and are very stable. That is, R is an appropriate quantity for 
the determination of 03, 

There are other theoretical uncertainties from higher-order 0(a;^) and higher-twist 0{1/Mb) 
corrections. For a simple estimation, we examine the fractional contribution to the form factor 
F^'^ as a function of as{t)/7T. It is observed that 90% and 97% of the contributions arise 
from the region with as{t)/7c < 0.2 and with as{t)/7[ < 0.3, respectively. Therefore, our 
PQCD results are well within the perturbative region. It is reasonable to assume that 0{al) 
corrections to the decay amplitudes are below 20%. In the derivation of the hard functions, 
we have neglected the mass difference A = Mb — mb to obtain the leading-twist factorization 
formulas. Next-to-leading-twist corrections, proportional to A/Mb, are then about 10%. 



5 PENGUIN ENHANCEMENT 

In this section we shall highlight the enhancement of penguin contributions observed in the 
PQCD approach, and its role in the explanation of the B ^ Tin and B Kir data. For 
simplicity, we demonstrate our observation by means of the FA approach. Consider the ratios 
R in Eq. (|^) and Rt^ defined by 

Br(i?g ^ K^TT^) 
Br(50 7r±7rT) 



= p^/ . ' (78) 



which can be written as 



aj^ + 2a xX Rb cos 4>3 

R = 2 ' I'^J 
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_ ax + 2aKX^Rb cos (j)3 . . 

" " Xm,[R, + 2a^{R,-cos<Ps)] ' ^ ' 

The factors 

ax = , = , (81) 

«! ai 

being negative values, represent the ratios of the penguin contribution to the tree contribution 
in the Ktc and tctc modes, respectively. It is obvious that the data i? ~ 1 imply ^3 ~ 90°, no 
matter what ax, A and Rb are. It is the reason when we vary all the parameters in the analysis 
in Sec. IV, the extraction of 03 remains invariant. 

While to determine 03 from the data of the ratio Rj^ ~ 4, one must have precise information 
of ax and a^r, and of the parameters A and Rb. It can be shown that the extraction of 03 
from depends on these parameters sensitively. Hence, Rt^ is not an appropriate quantity for 
the determination of 03. To explain the data of -R^r ~ 4 in the FA approach, an unreasonably 
large mo ~ 4 GeV corresponding to = 2m„ = 3 MeV, i.e., large |ai^(,r)| ~ 0.09 and a large 
03 ~ 130° must be postulated |jT9|. This is obvious from Eq. (|80D , since a large |a;^(7r)| enhances 
i?7r; and a large 03 leads to constructive interference between the two terms in the numerator 
of i?7r- The determination 03 ~ 114° from global fits to charmless B meson decays |TP|, located 
between the two extreme cases 90° and 130°, is then understood. The result of 03 will become 
even larger, if reasonable quark masses = 2mu ~ 10 MeV are employed. The huge difference 
between 90° and 130° extracted from different data renders the determination of 03 in the FA 
approach less convincing. In the modified FA approach with effective number of colors N^^, a 
large unitarity angle 03 ~ 105° is also concluded . 

An interesting question is as follows. If we give higher weight to the extraction of 03 from 
R, which is more model-independent than that from R^^, can we explain the data -R^r ~ 4 using 
a smaller 03? The answer is positive in the PQCD approach. Table I shows that the ratio of 
the penguin contribution to the tree contribution reaches 

\ax\ = ^ - 0.1 , (82) 



even with a reasonable value of rriQ = 1.53 GeV. The reason is that we do not assume the same 
form factors for the operators Oi, 2,3,4 and for O^^q. These form factors, evaluated explicitly in 
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the PQCD formalism, possess different factorization formulas as shown in Eqs. ( [27| ) and (P^). 
It is easy to observe that the integrands in the two factorization formulas become identical, if 
the terms associated with the pseudoscalar wave function 0^ and the factors X3 are dropped. 
The X3 — limit corresponds to the kinematic configuration, in which the light quark emitted 
from the b quark decay vertex carries the full meson momentum. This is the configuration, 
on which the equality of the two form factors in the FA approach is constructed. Therefore, 
the larger ratio of the penguin contribution to the tree contribution is achieved dynamically, 
instead of by increasing mg. With this penguin enhancement, the observed branching ratios of 
the B — ^ Ktt and B ^ titt decays and i?^ ~ 4 can be explained simultaneously in the PQCD 
approach using mg = 1.53 GeV and a smaller 03 = 90°. That is, the data of i?^ do not demand 
large ttiq and 03. Such a dynamical enhancement of penguin contributions can not be obtained 
in the FA appraoch. 

One of the sources responsible for the penguin enhancement is the RG evolution effect caused 
by the running hard scale t. In Fig. 8 we display the RG evolution of the Wilson coefficients 
ai{^), i = 1,4,6. It is found that ai is almost constant for fi = 500 MeV to Mb- In contrast, 
1 04 1 and lag] dramatically increase as /i evolves to below Mb/2. If choosing t = Mb/2 with 
mo = 1.5 GeV, the ratio |ai^(7r)| ~ 0.06, close to that in the FA approach with the same value of 
mo, is too small to explain i?^ ~ 4. As stated before, PQCD provides a prescription for choices 
of the hard scales t: t should be chosen as the virtualities of internal particles in Eq. ( P7| ) in 
order to decrease higher-order corrections. It reflects the fact that energy releases and evolution 
effects involved in different B meson decay modes are different. These hard scales can then 
reach lower values, at which |a6(^)| is enhanced over \aQ{MB/2)\. This evolution effect increases 
^/^(Tr)! by about 50% as indicated by Eq. (^21). 

The enhancement due to the increase of 0^(1) with decreasing t makes us worry that the 
contribution from the small t region may be important. This will invalidate the perturbative 
expansion of the hard amplitudes. As a check, we examine the fractional contribution to as 
a function of as(t)/n. The results are displayed in Fig. 9, which indicate that about 80% (90%) 
of the contributions come from the region with as{t)/TT < 0.2 (0.3). Therefore, exchanged 
gluons are still hard enough to guarantee the applicability of PQCD. 
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Another source of penguin enhancement is the behavior of the B meson wave function at 



X ^ 0. As shown in Eqs. (27) and (|28|) , the factorization formulas consist of two terms. It 
can be easily verified that when the two terms are roughly equal, the ratio of the penguin 
contribution to the tree contribution reaches its maximum. A simple investigation reveals the 
approximate expressions of the hard functions at small momentum fractions, 

he{xi,X3,bi,b3) ~ ln(xiX3)lnx3 , 

he{x3, Xi, 63, 61) ~ ln(xiX3) Inxi . (83) 



A B meson wave function with other behaviors, say, 0^ ~ x or ^/x as x — > 0, leads to 
the dominance of the second term, and the penguin contribution becomes relatively smaller. 
While the B meson wave function in Eq. (0), which vanishes like as x — > 0, renders the 
contributions from the above two terms approximately the same. The penguin contribution 



corresponding to Eq. (|59|) is about 10% larger than that corresponding to the model in p? 



6 FINAL STATE INTERACTION 



Two-body FSl effects have been studied in various ways It was found that these effects 
enhance the CP asymmetry in the B^ K^ir^ modes from order 0.5 % under the FA |^ up 
to order (10 - 20) %. However, Kamal has pointed out that the large CP asjnnmetry is due to 
an overestimation of FSI effects by a factor of 20 [H6|. For a critic assessment on the analyses 



of FSI effects in the literature, refer to |^| 



We briefly sketch the methods used in most of the estimates of FSI effects. For simplicity, we 
consider only the B~^ — > K^tt~^ decay. The unitarity relation for the amplitude A{B^ K^n^) 
is written as 

Im^(S+ ^ K^7i+) = -Y. '^T^K^B - En)A{N K%+)A*{B+ N) . (84) 

2 N 

If only the elastic channel K^tt^ contributes, Watson's theorem tells that the phase of A{B^ 
K^TT^) is given by the S wave / = 3/2 phase shift. This argument works for the K tttt 
decays but not for B meson decays. For Mb ~ 5 GeV, many channels contribute and Watson's 
theorem says nothing about the strong phase of A{B^ — ^ K^n^). In fact, even if the phases 
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of ^(A^ — > K^TT^) for all are known, the unitarity relation does still not fix the phase of 
A{B+ K^Ti+) uniquely. 

Inspite of this difficulty, some authors computed ^(A^ — > K^tt^) for few A^. Certainly, more 
than the unitarity relation is needed to obatin the strong phase of A{B^ K^tt^). The phases 
of ^(A^ —>■ K^TT^) are often estimated by a Regge analysis. However, this method is reliable 
only near the forward direction. In our problem we need S wave amplitudes, i.e., scattering 
amplitudes for all angles. A big assumption of a straight line trajectory has been adopted. This 
is highly questionable, especially for Pomerons. For these reasons, we believe that the above 
analyses are qualitative at most. 

It is our viewpoint that if a strong phase can not be determined in QCD, there is no other 
way to compute it. A simple physical picture of FSI, the color-transparency argument, has 
been put forward by Bjorken ||47|| : 

Since products of a -B meson decay into two light mesons are quite energetic, 
the quark-antiquark pair inside a meson remains a state of small size with a corre- 
spondingly small chromomagnetic moment until it is far from the other meson. It 
is then more realistic that the two quark pairs group individually into final-state 
mesons without further exchanging soft gluons. 

This picture is consistent with our observation: Sudakov suppression is strong for large meson 
momenta as shown in Eq. (JT^), which then demands final-state mesons of small transverse 
extent b, i.e., of small chromomagnetic moment. 

The effects from soft gluon exchanges among mesons in two-body heavy meson decays 
have been analyzed quantitatively by means of RG methods, which sum up large logarithms 
produced by infinite gluon emissions. It was found that these effects generate only small FSI 
phases for B meson decays, in agreement with the color-transparence argument, but large FSI 
phases for D meson decays . That soft gluon effects are large in D meson decays is expected, 
since Sudakov suppression is weaker, two quarks in a final-state meson is separated by a larger 
distance, and soft gluons can resolve the color structure of that meson. Based on the above 
reasonings, we have neglected FSI effects in the PQCD approach to two-body B meson decays. 

To justify the neglect of FSI, we apply our formalism to the B ^ nn and B Kn 
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decays without taking into account these long-distance effects. FSI in these decays should be 
different. Since tree contributions dominate in the B tttt decays, extra phases from FSI do 
not change branching ratios very much. This argument applies to the decays — > K^tt^ 
and B^ — > K^tt^, where penguin contributions dominate. While the inteference between tree 
and penguin contributions plays an essential role in the B^ — > K^tc^ and B^ — > K^tt^ decays. 
Large FSI effects will change the relative phases between tree and penguin amplitudes, and 
thus branching ratios. If the same formalism without including FSI can be applied to both 
decays successfully, we believe that these long-distance effects are negligible. The agreement of 
our results with the data shown in Eq. ( |7l| ) implies this conclusion. 

Without FSI, large strong phases come from the factorizable annihilation diagrams (phases 
from nonfactorizable diagrams are small) in the PQCD approach as stated in Sec. IV. There 
has been a widely spread folklore that the annihilation diagrams give negligible contribution 
due to helicity suppression, the same as in vr — > eu decay. That is, a left-handed massless 
electron and a right-handed antineutrino can not fly away back to back because of angular 
momentum conservation. However, this argument does not apply to F^. A left-handed quark 
and a left-handed antiquark, for which helicities are dictated by the Oq operator, can indeed 
fly away back to back |4^. These behaviors have been reflected by Eqs. (|29|) and (pOl) : Eq. (^) 
vanishes exactly, if the kaon and pion wave functions are identical, while the two terms in 
Eq. ( pOD are constructive. The reason the annihilation diagrams from the Og operator possess 
large absorptive parts can be understood in the following way. The cuts on the internal quark 
lines in Figs. 3(e) and 3(f) correspond to a process 

B+ ^su^ K%+ . (85) 

The intermediate state (su) can be regarded being highly inelastic, if expanded in terms of 
hadron states. According to Eq. (|8^) , large strong phases are expected. 

7 COMPARISION WITH OTHER ANALYSES 



Recently, Beneke et al. proposed to evaluate nonfactorizable contributions to charmless B 
meson decays in the PQCD framework ^0|. They argued that factorizable contributions (tran- 
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sition form factors) are not calculable in perturbation theory, but nonfactorizable contributions 
are. The reasoning is as follows. The internal b quark in the hard amplitude may go onto mass 
shell, producing a power divergent factor x being the momentum fraction associated with 
the pion. The soft divergence from x ^ can not be removed by a pion wave funciton, if it 
vanishes like x as a; 0. Since this divergence is not of the pinched type which is absorbed 
into a wave function, its appearence implies the breakdown of PQCD factorization theorem. 



While such a power divergence does not exist in nonfactorizable amplitudes |5^. We argued 
that the factor in fact can be easily smeared out by parton transverse momenta kx con- 
sidered in this work or killed by a wave function vanishing faster than x as x ^ 0. With the 
inclusion of kx and Sudakov suppression, we have explicitly shown that almost 100% of the 
full contribution to the B ^ tt transition form factor arises from the region with the coupling 
constant a^/vr < 0.3. It indicates that dynamics from hard gluon exchanges indeed dominate 
in the PQCD calculation. 



There are other important differences between our approach and [^. The momentum of the 
light spectator quark in the B meson has been ignored in the formalism of |5^, such that quark 
propagators in hard nonfactorizable decay amplitudes always remain time-like. The annihilation 
diagrams were not included either. With these approximations, leading-order information of 
strong phases was lost. Strong phases then arise from diagrams of the BSS mechanism, which, 
as shown below, are small compared to those from annihilation diagrams. On the other hand, 
Sudakov resummation of large logarithmic corrections was not taken into account. It is then 
expected that higher-order corrections will be large and spoil the perturbative expansion. It 



has been shown |^ that the PQCD formalism without including Sudakov suppression is not 
applicable to exclusive processes for energy scale below 10 GeV. 

We show that strong phases from the BSS mechanism are suppressed by the charm mass 
threshold and by 0{as), since there must be a hard gluon emitted by the spectator as shown 
in Fig. 10, which turns the soft spectator in the B meson into a fast spectator in the final-state 
meson. That is, the contributions from the penguin diagram have been overestimated. The 
charm quark loop contributes an imaginary part, 

C2{t)as{t) J duu{l - u)e{q^u{l - u) - ml) , (86) 
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where is the invariant mass of the gluon emitted from the penguin. The contribution from 
the u quark loop is suppressed by the small CKM factor \Vu\- Since is not clearly defined in 
the FA approach, it is usually chosen as = ml or = ml/2, and Eq. ( p6D gives a substantial 
amount of imaginary contribution to decay amplitudes . 



However, the invariant mass q"^ can be defined unambigiously in the PQCD formalism by 

= {X2P2 + X3P3)' = X2XsMl , (87) 

since the quark going into the kaon (pion) carries the fractional momentum 0:2^2 [xsPs)- Then, 
= ml or ml/2 corresponds to a configuration, in which the two quarks produced from the 
gluon carry the full momenta of the two final-state mesons. Obviously, this configuration is 
unlikely because of the strong suppression from the kaon and pion wave functions in the large 
X region. Substituting Eq. ( pTl) into Eq. (^), an exact numerical analysis indicates that the 
BSS mechanism contributes an imaginary part smaller than that from the nonfactorizable and 
annihilation amplitudes by a factor of 10. Table II shows how the imaginay part of the charm 
quark loop contribution vanishes with the decarease of g^. 

On the issue of FSI, Suzuki has argued that strong phases of the B Ktc amplitudes 
can not be evaluated in QCD [Q. He pointed out that the invariant masses of the su and 
du pairs in Figs. 3(c) and 3(d), respectively, are of order (Aqcd^b)^^^ ~ 1-2 GeV. It implies 
that the B —* Kit decays are located in the resonance region and their strong phases are very 
complicated. We have computed the average hard scales of the B — > Kvr decays, which are 
about 1.4 GeV, in agreement with the above estimate. However, the outging quark pairs possess 
an invariant mass larger than 1.4 GeV, such that the processes are in fact not so close to the 
resonance region. We could interpret that the decays occur via a six-fermion operator within 
space smaller than (1/1.4) GeV~^. While they are not completely short-distance, the fact that 
over 90% of contributions come from the x-h phase space with as{t)/TT < 0.3 allows us to 
estimate the decay amplitudes reliably. We believe that the strong phases can be computed up 
to about 20% uncertainties, which result in 30% errors in the predictions for CP asymmetries. 
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8 CONCLUSION 



In this paper we have analyzed the B Kn decays using PQCD factorization theorem. In 
this approach hadronic matrix elements, including factorizable and nonfactorizable, and real 
and imaginary contributions can be evaluated explicitly. The strong phases arise from non- 
pinched singularities of quark and gluon propagators in annihilation and nonfactorizable hard 
amplitudes. It has been explicitly shown that strong phases from the BSS mechanism are small. 
The analysis of soft gluon effects and the simultaneous success of the PQCD applications to 
the B — >■ Kt: and B tttt decays implied that long-distance FSI effects are negligible. The 
universal meson wave functions have been determined from the available data of the pion form 
factor and of the B — > Dn and S — > tttt decays. The dependences of the ratio R of the neutral 
B decay branching ratio to the charged B decay branching ratio and of the CP asymmetries 
on 03 have been derived. Our predictions for all the B — > Kit modes are consistent with the 
experimental data. 

In spite of potential theoretical uncertainties, we have extracted the following features for 
the B — > Ktt, tttt decays, which are less ambiguous: 

1. Nonfactorizable amplitudes are negligible. 

2. Annihilation diagrams are not negligible. 

3. Annihilation diagrams generate large strong phases. 

4. More precise data are needed in order to obtain a strong constraint on 03 . 

5. is an ideal quantity for the determination of 03, since it is insensitive to all the Wolf stein 
and nonperturbative QCD parameters. 

6. 03 is about 90" from fitting our predictions to the central value of the data of R. 

7. Penguin amplitudes are dynamically enhanced, and larger than those employed in the FA 
approach by 50%. 
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8. The data oi B ^ rni decays, ie., the ratio R.,^ of the B Ktt branching ratio to the 
S — > TTTT branching ratio can be explained by the smaller angle 03 ~ 90°. That is, the 
data of do not demand a large 03 > 90°. 
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APPENDIX A 



In this Appendix we supply the details of the Wilson evolution. The matching conditions 
at /i = Mw are given by 



C2{Mw 
Cs{Mw 



0, 
1, 



-Eo{Xt) + 



247r 
as{Mw) 
Stt 



6n sin^ Q\y 



-[2Bo{xt) + Coixt)], 



247r 



Eoixt), 



6T[ 



OLf. 



0, 



\^C^{xt) \ D^{xt)\ 



Uc^{xt) + D^{x:t) + - 

I SI 



1 



sin^e 



w 



[105o(xt) - 4Co(a;t)] 



0, 



with Xt = m^/M^, rrit being the top quark mass. The functions Bq,Co,Dq and Eq are the 



Inami-Lim functions ||53 

Bo{x) = 
Coix) = 
Do{x) = 
Eo{x) = 



X x\nx 
1 — X (x — 1)^ 
X — 6 3x + 2 



X — 1 (x — 1)^ 



In X 



4, 19x^-25x2 x2(5x2-2x-6), 

In X ; — 1 ; -; lu X , 

9 36(x-l)3 18(x-l)4 

2, x2(15- 16x + 4x2) , x(18-llx-x2l 

In X ; -; In X H ; — 

3 6(x-l)4 12(1 -x)3 



(89) 
(90) 
(91) 
(92) 



We adopt the following parameters: mt = 170 GeV, Mw = 80.2 GeV, a^(Mvy) = 0.118, 



= 1/129, sin^ Qw = 0.23 and A^^ = 
The solution to Eq. (|58D is written as 



250 MeV. 



(93) 
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The evolution function including electroweak penguins is 

'-'em / 

Jin mi, 



U{t,mb,ae 



Uf{t,mh) 



rfln/x't/;(t,/x')[7f^];t/^(/x' ,'mb 



OLe. 



with 



Uf{t, rrib) = exp 



rhlt 

/ din// 

J In mi 



r-(0)Ti 
]/ 



For n = mb~ 4.8 GeV, the values of Ci{mb) are 



Ci(mb) = -0.271, = 1.124 , 

C3(mb) = 1.255 X 10"^ , C^imb) = -2.686 x lO'^ , 

C5{mb) = 7.805 X 10"^ , = -3.287 x 10"^ , 

Crimb) = 3.453 x 10"^ , Cs(mb) = 3.177 x 10"^ , 

Cg{mb) = -9.765 x 10"^ , Cioirrib) = 2.240 x 10"^ . 



(94) 



(95) 



(96) 



Values of the Wilson coefficients at different energy scales (i ~ 1.0 GeV, 1.5 GeV, 2.0 GeV, 2.5 
GeV, 3.0 GeV and 4.8 GeV are listed in Table III. 

APPENDIX B 



Below we summarize the parameters we have adopted in the numerical analysis of this work: 



1. Masses, decay constants and lifetimes: 

M^± = 0.139 GeV, 
Mk = 0.49 GeV, 
mu^4:.5 MeV, 
THs = 100 MeV, 
rub = 4.8 GeV, 
Mw = 80.2 GeV, 
U = 130 MeV, 
Tbo — 1.55 ps. 



M^o = 0.135 GeV, 
Mb = 5.28 GeV, 
md = 8.1 MeV, 
rric = 1.5 GeV, 
mt = 170 GeV, 
fs = 190 MeV, 
fx = 160 MeV, 
tb- — 1.65 ps. 
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2. QCD and electroweak parameters: 



Gf = 1.16639 X 10-5 GeV-2, 
as{Mz) = 0.117, 
A = 0.2196, 



A^^ = 250 MeV, 



MS 

aem = 1/129, 
A = 0.819, 







2 , ,2 A,2' 




\ UJb / 


1 2 



i?,, = Vp' + = 0.38. 

3. Meson wave functions: 

(I>b{x) — Nbx'^{1 — x)^ exp 

Nb 

(f)^{x) 
(f)'j^{x) 



91.7835 GeV , ujb = 0.4 GeV , 

^ :Ax(l -x)[l + 0.8(5(1 - 2x)2 - 1)] , 



V2iv: 

3 



3 
3 

%/2iVe 



f^x{l - x) , 

/kx(1 -x)[1 + 0.51(1 - 2x) + 0.3(5(1 - 2xf - 1)] , 
fKx{l - x) . 
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Figure Captions 



1. Fig. 1: Unitarity triangle and the definition of the angles 0j. 

2. Fig. 2: Feynman diagrams for the K'^n^ decays. 

3. Fig. 3: Feynman diagrams for the — > K'^tt^ decays. 

4. Fig. 4: Factorization of the B — > Kn decays in the PQCD approach. 

5. Fig. 5: Dependence of the branching ratios of the B — > Kn decays on 03 with the upper 
(lower) dashed fine corresponding to the B (B) meson decays. 

6. Fig. 6: Dependence of the ratio it! on 03. The dashed (dotted) fines correspond to the 
bounds (central value) of the data. 

7. Fig. 7: Dependence of CP asymmetries on 03. The dashed (dotted) fines correspond 
to the bounds (central value) of the data of the B^ — > K^n^ decays in (a) and the 
5° ^ K^TT^ in (b). 

8. Fig. 8: RG evolution of the Wilson coefficients aj(//), i — 1 ,A ,6, normalized by tfieir 
values at fj. = nib. 

9. Fig. 9: Frcation contribution to as a function of as{t)/7r. 
10. Fig. 10: Feynman diagram for an induced c (li)-quark loop. 
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5.577 X 10-^ 


Ff 


-5.537 X 10-2 


Fa' 


3.333 X 10-3 + i 3.181 x 10-^ 


Me 


-0.942 X 10-=^ + i 3.385 x 10"^ 


Mf 


2.931 X 10-^ - i 1.304 x 10"^ 




-9.397 X 10-5 - i 1.918 x 10"^ 



Table 1: Contribution to the K~^7r decay from each form factor and nonfactorizable 

amphtude. 
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Re[G] 


Im[G] 


ml 


-0.760 


2.025 


ml/2 


0.139 


1.775 


m|/3 


0.912 


1.178 




1.288 


0.392 


ml /A 


1.162 


0.000 



Table 2: Real and imaginary parts of the charm quark loop contribution G{q^) — —4 / duu{l — 
u) \n\m1 — q^u{l — u))] in the BSS mechanism. 
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1.0 GeV 


1.5 GeV 


2.0 GeV 


2.5 GeV 


3.0 GeV 


4.8 GeV 




0.5139 


0.4208 


0.3626 


0.3275 


0.3034 


0.2552 




-0.650 


-0.510 


-0.435 


-0.385 


-0.349 


-0.271 


C2 


1.362 


1.268 


1.219 


1.189 


1.168 


1.124 


Cs 


0.036 


0.027 


0.022 


0.019 


0.017 


0.013 


C4 


-0.063 


-0.050 


-0.043 


-0.038 


-0.035 


-0.027 


C5 


0.015 


0.013 


0.012 


0.011 


0.010 


0.008 


Ce 


-0.102 


-0.074 


-0.060 


-0.051 


-0.045 


-0.033 




0.040 


0.035 


0.035 


0.036 


0.038 


0.045 


Cs/ 01 em 


0.128 


0.091 


0.073 


0.062 


0.055 


0.041 


C'g/ CKg^ 


-1.509 


-1.416 


-1.366 


-1.334 


-1.311 


-1.260 


C'lo/ Oiem 


0.695 


0.546 


0.465 


0.412 


0.373 


0.289 



Table 3: Values of the running coupling constant as and the Wilson coefficients Ci with A 
250 MeV for different energy scales // = 1.0, 1.5, 2.0, 2.5, 3.0, and 4.8 GeV. 
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